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We study, using a perturbative renormalization group technique, the phase diagrams of bond- 
aligned and diagonal Hubbard ladders defined as sections of a square lattice with nearest-neighbor 
and next-nearest-neighbor hopping. We find that for not too large hole doping and small next- 
nearest-neighbor hopping the bond-aligned systems exhibit a fully spin-gapped phase while the 
diagonal systems remain gapless. Increasing the next-nearest-neighbor hopping typically leads to 
a decrease of the gap in the bond-aligned ladders, and to a transition into a gapped phase in 
the diagonal ladders. Embedding the ladders in an antiferromagnetic environment can lead to a 
reduction in the extent of the gapped phases. These findings may suggest a relation between the 
orientation of hole-rich stripes and superconductivity as observed in I^-^Sr^CuCU . 



c3 



i 

C 

O 
o 



(N 
> 

m 
o 
o 

o 

03 



o 
o 



X 



I. INTRODUCTION 



There is a growing number of experiments which 
suggest that the existence of inhomogeneous electronic 
structures is a generic feature of the hole-doped high- 
temperature superconductors. Often these structures 
take the shape of quasi-one-dimensional hole-rich stripes 
running along various directions in the copper-oxygen 
planes 1 . The relevance of inhomogeneities in general 
and stripes in particular to the phenomenon of high- 
temperature superconductivity is presently far from be- 
ing a settled issue. 

Neutron-scattering experiments have produced clear 
evidence for a significant correlation between signatures 
of incommensurate spin excitations in L^-zSr^CuC^ 
and the superconducting properties of this compound. 
Recently, it has been demonstrated 2 that such low-energy 
incommensurate spin excitations exist in the supercon- 
ducting phase of L^-zSr^CuC^ all the way up to its 
overdoped boundary at x = 0.28, where they disap- 
pear simultaneously with superconductivity. On the 
underdoped side of the phase diagram static magnetic 
stripe signatures exist both in the insulating spin-glass 
phase 3-5 with 0.02 < x < 0.05 and in the underdoped 
superconductor 6-8 with 0.05 < x < 0.12. However, the 
transition from the insulator to the superconductor at 
x = 0.05 is accompanied by a rotation of the stripe sig- 
nal from the diagonal direction in the former to a bond- 
aligned signal in the latter. Moreover, angle resolved pho- 
toemission spectroscopy (ARPES) has revealed 9 an ap- 
parent coincidence between the collapse of the gap along 
the nodal (diagonal) directions in A:-space and the ap- 
pearance of diagonal stripes in samples near the Neel 
boundary of the antiferromagnet at x = 0.02. 

Motivated by these observations we are interested in 
finding a possible correspondence between the geome- 
try of stripes and their low-energy properties. However, 
studying the physics of stripes as self-organized systems 



in the strongly interacting copper-oxygen planes is a diffi- 
cult task for which reliable analytical tools are still miss- 
ing. Instead, we propose to investigate rudimentary mod- 
els of stripes in the form of Hubbard ladders defined 
as bond-aligned and diagonal sections of a square lat- 
tice with nearest-neighbor t and next-nearest-neighbor t' 
hopping. In doing so we avoid the question of the mech- 
anism which leads to the segregation of holes into stripes 
(on which there are different views 10-13 ), and assume 
that it takes place at energies higher than the energy 
scales which govern the physics of holes on the ladders. 
On the other hand, we will concentrate on a single lad- 
der thus restricting the discussion to energies larger than 
any inter-stripe coupling (such as inter-stripe tunneling) , 
which is assumed small. While the degree to which holes 
are restricted to move along stripes in the real systems 
is not well understood it is clear that the limit of total 
confinement, assumed in the ladder models, is an ideal- 
ization. To probe the sensitivity of the results to this 
point we consider, for each of the two geometries men- 
tioned above, both 2 and 3-leg ladders. 

Arguably, the most severe shortcoming of the present 
study, with respect to its applicability to the physics of 
stripes in the cuprates, is the fact that the interaction 
between electrons on the ladders is assumed weak and 
treated using a perturbative renormalization group (RG) 
technique. Nevertheless, previous numerical studies 14-19 
have shown that many of the qualitative features uncov- 
ered by perturbative analysis of Hubbard ladders sur- 
vive the limit of strong coupling where the emergent en- 
ergy scales are of the order of the exchange interaction J. 
This fact suggests that the gross features which we dis- 
cover in the phase diagrams of the models considered here 
may likewise extend into the strong coupling regime. Of 
course, this needs to be verified by explicit calculations. 

Our principal results are presented in Figs. 2-5. They 
contain the phase diagrams of the ladders as a function 
of the relative strength of next-nearest-neighbor hopping 
t'/t and the average electronic occupation per site n. 



1 



The various phases are labeled according to the num- 
ber of charge and spin modes 20-22 that remain gapless 
in the presence of interactions. We find a distinct differ- 
ence between the behavior of bond-aligned and diagonal 
ladders for small values of t'/t and doping levels in the 
range of0.5<n<1.5. In this region the bond- aligned 
ladders are maximally-gapped with only one remaining 
gapless charge mode associated with global gauge and 
translational invariance 22 (at commensurate fillings this 
mode may be gapped as well) . This finding extends pre- 
vious results for such models obtained both from weak- 
coupling treatments 20-26 and numerical calculations at 
strong-coupling 14-18 . Conversely, the diagonal ladders 
remain gapless under similar conditions. Such behavior 
has been known for the 2-leg "zig-zag" ladder 27 ' 19,28 and 
our results demonstrate that it persists also in the case 
of the diagonal 3-leg "diamond" ladder. 

Increasing the value of t'/t in bond-aligned systems, 
which exhibit a gapped phase for zero t' , leads to a de- 
crease of the gap. In the diagonal ladders similar vari- 
ation of the next-nearest-neighbor hopping results in a 
transition from a gapless to a gapped phase above a cer- 
tain critical value. Further increase of t'/t enhances the 
gap, which reaches a maximum and then decreases. 

In order to gain some insight into the effects of a mag- 
netically ordered background on the physics of ladders we 
replace the non-interacting spectrum of the ladder by the 
bound-states spectrum of a linear defect in an antiferro- 
magnetic two-dimensional lattice. In this way we include 
on a crude mean-field level the interactions between the 
stripe electrons and the surrounding spins. The inter- 
actions between electrons on the defect are still treated 
perturbatively. We find that the changes in the spec- 
trum typically have little effect on the overall structure 
of the phase diagrams, as long as the electronic states are 
reasonably localized on the defects. When changes do oc- 
cur they often tend to reduce the extent of the gapped 
phases. 

Notwithstanding the above mentioned limitations of 
our study, we note that our findings, and the experi- 
mental congruity between stripes orientation and super- 
conductivity, conform with the notion that the origin of 
the superconducting gap is to be found in the gap which 
develops on the quasi-one-dimensional stripes 24 . The re- 
sults may also provide an explanation, within the stripes 
picture, for the dichotomy between gapped and gapless 
behavior in fc-space, as revealed by ARPES. We elabo- 
rate on these points and raise few difficulties related to 
such interpretations in the discussion section. 



II. STRIPES AS LADDERS 

A. Models and methods 

Consider the Hubbard model on a square lattice with 
nearest-neighbor t and next-nearest-neighbor t' hopping. 



We begin by modeling a single stripe as a quasi-one- 
dimensional section consisting of N connected infinite 
chains, i.e. an N-leg ladder, on this lattice. In this pa- 
per we will concentrate on the minimal models for bond- 
centered and site-centered stripes as 2-leg and 3-leg lad- 
ders, respectively. We will consider two geometries for 
the ladders: bond-aligned and diagonal. In the first the 
ladder lies along the nearest-neighbor bonds while in the 
second it extends in the direction of the next-nearest- 
neighbor bonds (See Fig. 1). In both cases we assume 
open boundary conditions in the transverse direction. 
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FIG. 1. 2 and 3-leg ladders on the background lattice, (a) 
bond-aligned ladders (b) diagonal ladders. 

The bond-aligned N-leg ladder is described by the 
Hamiltonian H = H% ond + H v 



N 



H, 



bond 



EE E 

i=\ x a=tl 
t[ C i,X,a C i+l,x,<j + C i,x,cr C i,x+l,cr + h.c] 

t l C i,x,cr C i+l,x-l,cr + C i,x,a C i+l,x+l,cr + h.C.] 
N 



(1) 

(2) 



Hu = ^EE 1 n i^A n i^,l : ' 
i— 1 x 

where c i x a annihilates a spin a electron at site x on 
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chain i and rii tX}(T = c\ ty c i x a . For the bond-aligned lad- 
ders we take as unit length the distance between nearest 
neighbors. The on-site Hubbard interaction is assumed 
weak [/<f. 

The diagonal N-leg ladder is described by the Hamil- 
tonian H = H^ ag + H v where 



N r 

EEE { 

i=l x cr=T,4. k 



*[ c i,x,cr C i+ i ;X+ i a + c i,x,<7 C i + i iX _i a 



h.c] 



- *'[ C lx,a C i,x+l,a + C L- C *+2,*, CT + fe-C.]| • (3) 

For this model we choose the next-nearest-neighbor dis- 
tance as our unit of length along the x direction. This 
means that x runs over the integers on odd chains and 
over the integers shifted by 1/2 on even chains. 

We study the above models in the weak-coupling limit 
using RG techniques. Much of the analysis is identical to 
the one used in previous studies of Hubbard ladders 20-26 
and consequently we outline here only the basic steps of 
the procedure, relegating some of the remaining details 
to the appendix. 

Since the interaction term is treated as a perturbation 
the first step is the diagonalization of the quadratic part 
of the Hamiltonian via the transformation 



N 

E z. 

b—l k x — — i 



1- h. j- .l 



Sibi>b,a(k x ) 



(4) 



where the matrix S depends on the specific model. The 
resulting free band structure is then rilled up to the chem- 
ical potential (we assume zero temperature throughout 
the paper) which is determined by the averaged electronic 
occupation per site n. Each band that crosses the chem- 
ical potential gives rise to Nf, symmetric pairs of Fermi 
points, where more than one pair may exist for a sin- 
gle band, see Fig. 8 for an example. A linearized chiral 
electronic branch is associated with every Fermi point. 
In order to have positive definite Fermi velocities we as- 
sign left (right) moving fields to Fermi points at which 
the spectrum is a descending (ascending) function of k x . 
When the expansion of the band field operators in terms 
of these chiral fields 

N b 

(5) 

with r\ = sign(fcp) at the right moving point, is 
plugged into the Hubbard term, Eq. (2), a host of 
interaction terms emerges. Typically only forward- 
and Cooper-scattering are allowed by momentum and 
spin conservation 20,21 , but for appropriate values of the 
parameters other momentum-conserving "special" pro- 
cesses and Umklapp terms appear as well (See Figs. 2,3 
and the appendix). 



In the first stage of the renormalization process the RG 
equations for the various coupling constants (as given in 
the appendix) are integrated numerically until I, the loga- 
rithm of the cut-off in units of its initial value, approaches 
a point I* where one or more of the couplings grows as 
(I — I*) i with 7 < 1. At this scale the flow has reached 
the vicinity of a strong-coupling fixed point. The most 
divergent couplings in this region tend to freeze certain 
modes and open a gap in their fluctuation spectrum. The 
pinned modes may be identified via bosonization and the 
gap scale associated with them may be related to the 
original bandwidth Ep according to 



Ai w E F e~ 



(6) 



To obtain the behavior of the system on scales larger 
than I* one needs to carry out a second RG stage in which 
the relevance of the weakly or non-diverging couplings of 
the first stage is assessed near the new strong-coupling 
fixed point 22 ' 26 . This is done by calculating their scaling- 
dimcnsion 5 while freezing the gapped modes. In case 
some of these couplings become relevant the flow is di- 
rected towards a new fixed point in which additional 
modes are gapped with a typical scale of 



1/(2-5) 



(7) 



where <?* is the dimensionless strength of the residual in- 
teractions at the end of the first RG stage 22 ' 26 . In prin- 
ciple one should repeat the process until a stable fixed 
point is reached. In our case we find that if new rele- 
vant couplings do occur in the second stage, they take 
the system into the maximally allowed gapped phase. 



B. Results 

Our results for the ladder models are summarized in 
the phase diagrams as presented in Figs. 2-5. To facili- 
tate comparison with the cuprates we concentrate on the 
case t' < 30 ' 31 . Since the model is defined on a bipartite 
lattice the spectrum is invariant under t — > —t. Moreover, 



the particle-hole transformation c x , y , 



maps the 

original Hamiltonian on itself with t — > —t, t' — > —t' 
and changes the average electronic density according to 
n —* 2 — n. One can therefore obtain the t' > half of the 
diagrams by reflecting them with respect to the t' = 
and the n = 1 lines. 

We use the CraSn classification scheme 20-22 wherein a 
phase is characterized by the number (m, n) of its gap- 
less charge and spin modes. In the diagrams we shaded 
the regions in which the systems end up in the maxi- 
mally gapped phase (C1S0) allowed by global gauge and 
translational invariance 22 . In the lightly shaded areas 
this phase is obtained only at the end of the second RG 
stage and the nature of the phase at the end of the first 
stage is shown in parentheses. More information on the 
character of the massive fields, and the relevant operators 
which give them their masses, is detailed in the appendix. 
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We find that whenever the CI SO phase is not reached the 
systems retain all of their non- interacting gapless modes, 
as indicated in the clear parts of the diagrams. 

n 



0.5 1 1.5 




C1S1 



FIG. 2. The phase diagram of the 2-leg bond-aligned lad- 
der. The system remains gapless in the light regions and ex- 
hibits the maximally gapped C1S0 phase in the shaded parts 
of the diagram. In the darker shaded regions this phase is 
obtained at the end of the first RG stage while it is reached 
only at the end of the second stage inside the lightly shaded 
domains. In the latter case we indicate in parentheses the 
nature of the phase at the end of the first stage. Umklapp 
and "special" processes are allowed on the dashed lines. The 
number of non-interacting Fermi points changes by 2 across 
the thick solid lines and the Fermi velocity associated with 
the pair vanishes there. On the dotted line a flat band oc- 
curs. Consequently, our results are unreliable in the vicinity 
of these lines. 

For the 2-leg ladders we have also considered the possi- 
bility of Umklapp scattering and momentum-conserving 
"special" processes which exist for particular values of 
the Fermi wave- vectors. The dashed lines in the diagrams 
correspond to parameters for which such processes are al- 
lowed. In all cases we find that the nature of the resulting 
phases, as indicated next to the lines, is determined by 
the end of the first RG stage. Along the sections marked 
by "?" the relevant terms of the first RG stage tend to 
pin pairs of canonically conjugated fields. Consequently, 
we were unable to determine the fate of the system in 
this eventuality. 

Our renormalization procedure relies on the lineariza- 
tion of the non-interacting spectrum around the Fermi 
points. It therefore fails whenever a pair of Fermi points 
resides at an extremum of the band structure. This hap- 



pens along the thick solid lines in the diagrams. Across 
each of these lines the number of Fermi point pairs 
changes by 1 (or 2 in the case of a "W" shaped band). 
For the bond-aligned ladders a flat band arises along the 
dotted lines in their diagrams while for the 3-leg diagonal 
ladder such a band occurs for t' = and 2/3 < n < 4/3. 
Our results are unreliable in the vicinity of these lines. 

Our most important finding is a correlation between 
the ladder orientation and the nature of its phase at 
small values of t' /t. The bond-aligned ladders exhibit 
the gapped CI SO phase over a wide range of electronic 
filling 0.5 < n < 1.5. The diagonal ladders on the other 
hand remain gapless, except for a small region near n = 1 
in the case of the diagonal 3-legged system. 

n 



0.5 1 1.5 




FIG. 3. The phase diagram of the 2-leg diagonal ladder. 



Bond-aligned ladders with t' = have been analyzed 
in the past using perturbative RG 20-26 and in most cases 
were found to be in the C1S0 phase 32 . This conclusion 
was further supported by density matrix renormaliza- 
tion group (DMRG) calculations at strong-coupling 14-18 . 
Our results demonstrate the robustness of this phase in 
the presence of not too large next-nearest-neighbor hop- 
ping. Our results also agree with the gapless behavior 
found in a previous DMRG study of the diagonal 2-leg 
ladder 19 and add information about this system at com- 
mensurate fillings. We are unaware of past treatments of 
the 3-leg diagonal ladder, especially in the strong inter- 
action regime. It is therefore unclear whether the gapless 
phase found by us in this model survives beyond the limit 
of weak-coupling. It is, however, interesting to note that 
at half filling and for t 1 = and U — ► oo, the diagonal 
3-leg ladder is equivalent to an alternating spin-1 and 
spin-1/2 chain, which is known to be gapless 33 . 
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FIG. 4. The phase diagram of the 3-leg bond-aligned lad- 
der. To distinguish between gapped phases which originate 
from a C2S2 phase, and those coming from a C3S3 phase, we 
mark the latter by an asterisk. 



C1S0 phases, regardless of the RG stage at which they 
become fully gapped. One exception is the 3-leg bond- 
aligned ladder at small t'/t where the gap increases con- 
siderably in the region around n — 1 due to a much 
smaller I* (which, in passing, we note scales as t/U). 
The gap is very small in the sliver phases. 
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FIG. 6. Gap scale in the 2-leg bond-aligned ladder. 
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FIG. 5. The phase diagram of the 3-leg diagonal ladder. 

One may estimate the gap scale in the gapped phases 
from Eqs. (6,7). Typically we find that the gap changes 
smoothly across boundaries between adjacent extended 




FIG. 7. Gap scale in the 2-leg diagonal ladder. 



The dependence of A on t' for the 2-leg ladders is pre- 
sented in Figs. 6,7. Increasing the magnitude of the 
next-nearest-neighbor hopping results in a decrease of 
the gap in the bond-aligned system. On the contrary, 
the diagonal ladder, which is gapless for t' = 0, develops 
a gap above a critical density-dependent value of t' . Fur- 
ther increase of the hopping amplitude causes the gap to 
grow, reach a maximum and then decrease. We observe 
a similar behavior also in the 3-leg ladders. 
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III. STRIPES AS DEFECTS IN AN 
ANTIFERROMAGNETIC BACKGROUND 



A. Model 



Since the cuprate high-temperature superconductors 
are born out of parent antiferromagnets one expects to 
find slow fluctuations of a collective field representing the 
local staggered magnetization. Stripes may be thought 
of as defects in this field along which holes tend to seg- 
regate. In the previous section we have assumed that 
the confinement of the holes to the stripes is extremely 
strong such that their coupling to the antifcrromagnetic 
environment may be neglected. Here we would like to 
relax this constraint. To do so we consider a model of 
non-interacting electrons on the square lattice interact- 
ing with a static staggered field 



defect \ ^ \ " 
cr=t,l 



x,y,(J x-\-X,y,a 



+ c. 



x ,y ,cr x ,y-\-l ,<r 

t 



h.c] 



y+l,a + c x,y,cr c x+l,y+l,a + h.C.] 



- r\r* r 
u V^x,y,a'-x — l 



(8) 



where c x y a annihilates a spin a electron at site (x, y) on 
the lattice (as defined in the coordinate system of Fig. 
la.) The field configuration is taken to represent defects 
with the same geometries as the ladders studied in Sec- 
tion II 



to left of ladder 
m(x, y) = I on the ladder 

±to right of ladder . 



(9) 



We consider both in-phase (+m) and anti-phase (— m) 
domain walls. 

In the absence of a defect, i.e. m(x,y) — to, 
the spectrum exhibits two bands of extended Bloch 
states whose energies E = — 4t' cos(k x ) cos(k y ) ± 
^/4i 2 [cos(fc I ) + cos(ky)] 2 + to 2 are separated by an en- 
ergy gap of order to for large enough magnetization. In- 
troducing a domain wall into the periodic background 
leads to the appearance of mid-gap states, see Fig. 8. 
These states are localized in the transverse direction to 
the wall 29 over a distance which typically scales as m -1 . 
We therefore expect their energies to approach, in the 
limit to — > oo , the non- interacting energies of the corre- 
sponding ladder. 
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FIG. 8. The band structure of the 2-leg bond-aligned de- 
fect for t'/t = —0.3. (a) The full spectrum including the 
bands of extended states for the case of an anti-phase defect 
with m/t — 5, (b) the localized bands for infinite m, (c) the 
localized bands of the anti-phase defect with m/t = 5, (d) the 
localized bands of the in-phase defect with m/t = 5. 

In the following we assume that m is large enough, 
such that varying the density of the system away from 
half filling is achieved by doping electrons or holes solely 
into the mid-gap bands. Consequently, we consider a 
reduced model whose non-interacting spectrum is com- 
posed of these localized mid-gap bands only. This band 
structure expresses, albeit on a coarse mean-field level, 
the effects of interactions between the electrons on the de- 
fect and the electrons comprising the antiferromagnetic 
background. In the cuprates these interactions are ex- 
pected to be of similar magnitude to the interactions be- 
tween the electrons moving along the stripes. However, in 
our perturbative treatment the latter are assumed much 
weaker. While unjustified when modeling real stripes 
this assumption enables us to treat the interactions on 
the defect in a controlled manner and gain some insight 
into the effects of the coupling between stripes and their 
antiferromagnetic environment. 

To make contact with the ladder models of the previ- 
ous section we relate the electronic operators on the de- 
fect sites to the creation operators of the localized states. 
The result is analogous to Eq. (4) with two differences. 
First, the exact expansion of the site operators gener- 
ally includes components along extended states. Such 
components are projected out and the remaining sum is 
properly normalized. As long as the mid-gap states are 
well localized on the defect the projected expansion re- 
mains a faithful representation of the site operators on 
the stripe. Secondly, for an iV-leg bond-aligned defect the 
staggered field induces a doubling of the unit-cell and of 
the number of localized bands. Accordingly, in Eq. (4), 
x becomes the unit-cell coordinate and i stands for the 
leg number together with the site index within the unit- 
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cell. The sum over the band index then runs from 1 to 
2N and the sum over the momentum along the defect ex- 
tends from — 7r/2 to it/2. This effect is absent in the case 
of diagonal stripes. On the other hand, and in contrast 
to the bond-aligned systems, the spin degeneracy of the 
bands is lifted whenever the spin texture is symmetric 
on opposite sides of a diagonal defect [this happens in 
anti-(in-) phase defects with even(odd) N}. 

From this point onward the analysis of the interactions 
between electrons on the defect, as given by the Hub- 
bard term Eq. (2), follows the same course described 
for the ladder models. However, owing to the absence 
of spin SU(2) symmetry in this case, we are compelled 
to consider, separately, scattering processes in different 
spin channels. 25 The resulting RG equations will be given 
elsewhere. 
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FIG. 9. The dependence of the phase diagram of an 
anti-phase 2-leg bond-aligned defect on the magnetization m. 
The results are for the case t' /t = —0.1. 



responding diagram for the analogous in-phase system, 
on the other hand, shows no significant changes in the 
phase boundaries down to the lowest value of m. A sim- 
ilar result is obtained for the 3-leg bond-aligned defects, 
with the exception that it is the in-phase system whose 
gapped phase is depleted at small m while the diagram of 
the anti-phase defect remains largely unaffected. We are 
therefore led to the conclusion that a larger sensitivity 
to changes in m exists when spins across the defect are 
aligned. 



-0.2 
-0.4 
-0.6 
-0.8 



-0.2 
-0.4 
-0.6 
-0.8 





n 




0.5 


1 


1.5 


(a) 




C1S1 


C1S0 


(C2Sl)-\ 
C2S2^\ 






n 




0.5 


1 


1.5 


(b) 










C1S1 


C1S0 


(C2S1) 

/ :iso \\ 

(C2S1) — 



FIG. 10. A comparison between the phase diagrams of (a) 
the 2-leg diagonal ladder and (b) the in-phase 2-leg diagonal 
defect with m — 3.5. 



B. Results 

In Fig. 9 we present the phase diagram of the anti- 
phase bond-aligned 2-leg defect as a function of m and n 
for a fixed small value of t' . We already referred to the 
fact that identifying the mid-gap sector of the spectrum 
with the electronic states of a stripe makes sense as long 
as these states are well localized on the defect. Since the 
localization length becomes larger with decreasing mag- 
netization or increasing next-nearcst-ncighbor hopping it 
imposes a limit on the minimal m and maximal t'/t which 
may be studied using this model. This constraint sets the 
lower bound of the magnetization axis in the figure. 

As can be seen lowering the modulation depth of the 
antiferromagnetic environment leads eventually to a re- 
duction in the extent of the gapped phase. The cor- 



The changes in the phase diagrams are initiated by de- 
formations in the non-interacting band-structure as m is 
reduced. For infinite m the mid-gap spectrum of the de- 
fect is identical to that of the corresponding ladder. For 
finite magnetization some of the degenerate points in the 
spectrum of the ladder separate by mini-gaps. When the 
background spin texture is symmetrical with respect to 
the defect the degeneracy at k x = ir/2 is lifted. When 
it is asymmetrical the bands separate at k x values away 
from the edge of the Brillouin zone, see Fig. 8. When a 
separation occurs, the result may be the removal of two 
Fermi points from the system (an example appears in 
Fig. 8c around E = —1). Such an event typically leads 
to a transition from a gapped to a gapless phase. In the 
symmetrical case, this scenario is more frequent, a fact 
which explains why it is more sensitive to variations in 
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to. Opposite occurrences where a band deforms from a 
"U" to a "W" shape, with the result of the appearance of 
an additional pair of Fermi points, are also possible but 
are less common. These changes lead to fragmentation of 
the phase diagrams of the defects. Usually gapped phases 
arc being penetrated by extensions of gapless phases, but 
we also encounter pockets or slivers of strong (first RG 
stage) C1S0 phase inside regions which are gapless or 
weakly gapped in the ladder system. 

We already mentioned that the spin degeneracy of the 
spectrum is lifted in the anti-phase 2-leg and the in- 
phase 3-leg diagonal defects. In these cases the RG anal- 
ysis indicates that the coupling constants do not flow 
and the systems remain in their gapless non-interacting 
phases. The spectra of the related systems with anti- 
symmetric spin configurations across the defects are still 
spin-degenerate and contain similar openings of mini- 
gaps as the bond-aligned defects. Consequently, their 
phase diagrams exhibit a more fragmented behavior, as 
demonstrated by Fig. 10. 

We find that lowering the magnitude of the staggered 
field diminishes the extent of the C2S2 phase throughout 
the phase diagrams. This effect is especially significant in 
the case of the 3-leg anti-phase diagonal defect where it 
leads to a reduction in the extent of the gapless phase in 
small t'/t, narrowing it down to a region centered around 
half filling. 



IV. DISCUSSION 

The results presented in this paper were derived for 
weak interaction strength U where all gap scales are ex- 
ponentially small. It is unclear to what extent they rep- 
resent the properties of ladders at strong coupling. How- 
ever, based on past experience coming from numerical 
studies 14 ^ 19 , it seems that many of the qualitative fea- 
tures of doped Hubbard ladders, such as the nature of 
their low-energy phases, are robust against changes in 
the interaction strength. In the following we assume that 
such continuity holds also for our results and, keeping in 
mind the other limitations of modeling hole-rich stripes 
as doped ladders, comment on their possible relevance to 
the physics of the cuprate high-temperature supercon- 
ductors. 

By now a large body of experimental data points to the 
fact that inhomogeneous electronic structures are preva- 
lent in the cuprates 1 . There is no consensus, however, 
on the role played by such inhomogencitics in the for- 
mation of high-temperature superconductivity and the 
other unique features of these compounds. One view, 
identifies the scattering of quasiparticles from the crit- 
ical charge fluctuations associated with a stripe quan- 
tum critical point, as responsible for the anomalous 
electronic properties, such as the pseudogap in the un- 
dcrdoped regime 34,35 . Another view, stipulates that 
the spin-gap, which develops on stripes, is the source 



of the pseudogap . The opening of this gap leads 
to intra-stripe power-law superconducting correlations, 
which turn eventually, via inter-stripe Joscphson pair- 
tunneling, into global long-range superconducting order. 

For stripes to play an essential role in the physics of the 
cuprates they must exist, either as static order or as fluc- 
tuating correlations, above the pseudogap crossover and 
the superconducting transition temperatures. Moreover, 
if the second point of view, just mentioned, is correct, 
these stripes must support a gap in order for supercon- 
ductivity to be achieved. 

As indicated in the introduction La2- a; Sr a ;Cu04 is a 
system where the existence of stripes is very well estab- 
lished over a wide range of hole-dopings and tempera- 
tures. As such it fulfills the first requirement. Our finding 
that under generic conditions weakly interacting diagonal 
ladders are gapless while bond-aligned ladders possess a 
spin-gap, offers some support to the hypothesis that a 
similar difference between the low-energy properties of 
bond-aligned and diagonal stripes is behind the coinci- 
dence of the insulator-superconductor transition and the 
rotation of the stripes direction in this compound 36 . It 
also supports the identification of bond-aligned stripes as 
the origin of the gapped ARPES spectrum around (ir, 0) 
(and symmetry related points) in fc-space, and the asso- 
ciation of the gapless Fermi surface at (7r/2,7r/2) with 
diagonal stripes. 

One may raise, however, difficulties associated with 
this picture. The first concerns signs of gapped com- 
ponents in the spin-glass phase of La2- ;t: Sr 2: Cu04 with 
x < 0.05. ARPES measurements have found a distinc- 
tive "knee" in the energy distribution curves at (tt, 0) in 
x = 0.05 samples 37-39 . This feature appears at a bind- 
ing energy of about 0.1 eV and becomes more prominent 
and less gapped with increased doping levels. It is con- 
sistent, within experimental uncertainties involving the 
doping, with the bond-aligned incommensurate signal ob- 
served in elastic neutron-scattering experiments 5 . This 
signal coexists with incommensurate diagonal peaks for 
0.056 < x < 0.07 and grows stronger with x. There is, 
however, some evidence coming from ARPES 37-39 that 
a similar but faint gapped band exists also in x = 0.03 
samples for which no corresponding neutron-scattering 
signal has been reported. 

A more serious objection, perhaps, may be raised in 
connection to the existence of a gapless diagonal compo- 
nent well inside the superconducting regime. A peak in 
the ARPES spectrum near (tt/2,it/2) and at the Fermi 
energy exists in all samples with 0.03 < x < 0.15 (and 
beyond) 39,38 . Its spectral weight intensifies upon enter- 
ing the superconducting region and continues to grow at 
least until optimal doping 39 . Although the above men- 
tioned coexistence of bond-aligned and diagonal incom- 
mensurate neutron-scattering peaks might explain this 
behavior near the insulator-superconductor transition, it 
does not offer an explanation at higher doping levels. One 
possible solution to the problem is that diagonal stripes 
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continue to exist deep in the superconducting phase, 
but only as fluctuating correlations. Such slowly fluc- 
tuating bond-aligned incommensurate correlations have 
been observed in inelastic neutron-scattering measure- 
ments of La 2 - a; Sr a; Cu04 with 0.12 < x < 0.27 2 ' 40 . 
Their detection along the diagonal directions would tie 
La 2 _a;Sr a ;Cu04 with another member of the cuprate 
family, YBa2Cu307_«5 , in which an anisotropic incom- 
mensurate ring was recently found by inelastic neutron- 
scattering experiments on undcrdoped and nearly opti- 
mally doped samples 41 . 
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APPENDIX A: THE RG PROCEDURE FOR THE 
LADDER MODELS 

1. Cooper- and forward-scattering 

As mentioned previously, spin and momentum conser- 
vation generically allows only for Cooper- and forward- 
scattering 21 . These processes are conveniently described 
using the charge and spin currents 

a 

^^E^^)^/^):, (Al) 

a,(i 

where er denotes the Pauli matrices. The same defini- 
tions apply for the left components. The allowed terms 
are 

J i,j=l L 
+ fij J ii J jj ~ fij J U ■ J jj^j i ( A2 ) 

where Vi is the Fermi velocity at ■ fij and Cjj are the 
forward- and Cooper-scattering amplitudes respectively. 
As fa and cn describe the same process we take fa = 0. 

The RG equations for these couplings and their bare 
values in terms of U and the transformation matrix S 
introduced in Eq. (4) have been derived by Lin, Balents 
and Fisher in Ref. 21. The most diverging couplings, as 
found by numerical integration of the RG equations, de- 
termine the nature of the phase on the scale at which the 
divergence occur. Subsequent scaling-dimension analysis 
determines the ultimate character of the system on even 
longer scales. The precise point at which a particular op- 
erator is declared most relevant is somewhat arbitrary. 



This uncertainty might affect the extent of the phases 
identified. Notwithstanding, we have verified the robust- 
ness of our results with respect to small variations in this 
criterion. Note that we chose not to display in the phase 
diagrams, Figs. 2-5, some of the phases whose ranges 
were found to be very small. 

When none of the couplings diverge the result is a gap- 
less phase. The scaling-dimension analysis indicates that 
all the couplings are marginal, thus ensuring the stability 
of the resulting CnSn phase (where n denotes the number 
of pairs of Fermi points). 

A CnSn-i phase is obtained at the end of the first 
rcnormalization stage when a single c a aa coupling di- 
verges. The bosonized relevant interaction is propor- 
tional to c" aa cos(v / 87T0a) an d results in a single gapped 
spin mode. A scaling-dimension analysis, after freez- 
ing (f>°, indicates that the interactions c°j cos[V2tt(9p — 
6j)] cos(v / 2~7r0j ) become relevant. They lead to a new 
C1S0 fixed point. 

The maximally allowed gapped phase, C1S0, is reached 
at the end of the first RG stage when all couplings except 
ffj diverge. When this happens c?(/*) w 4c? (/*) and the 
relevant interactions take the form 

xc& cos[\/2i(0£ - 9 p )] cos(v^tt#) cos(\/2^6J) . 

(A3) 

The terms with i = j pin all the spin modes. The terms 
with i 7^ j imply n — 1 gapped charge modes. The out- 
come is correspondingly a C1S0 phase. 

In the region where three pairs of Fermi points exist we 
encounter a situation in which all the couplings associ- 
ated with two Fermi point pairs (a, b) diverge at the end 
of the first RG stage, except for f% b . The relevant inter- 
actions are then similar to those of Eq. (A3), where i and 
j take the values a, b. The i = j terms lead to two gapped 
spin modes. The i j terms result in one gapped charge 
mode which is associated with 6 p ab ~ = (6? — 9 p )/\/2. 

The stability of the resulting C2S1 phase is deter- 
mined by the value of Kf£, the Luttinger exponent of 
Kb = ( 6 a + 0&)/\/2- When l/K p + < 4 the interactions 
associated with c a ac and c bc become relevant and the sys- 
tem ends up once again in the C1S0 phase. The value of 
this Luttinger exponent 20 

iKir 2 = (A4) 
4 ^ 2 - mi + y?r; b + ^ - i) 2 (Q 2 ] 

(1 + z,) 2 [2^ + {y - \yf; h \[2v - (1 + Qu + z, 2 )/^] ' 

where v — v a /vb, depends on f p £, the forward-scattering 
amplitude between the spin-gapped bands at the end of 
the first RG stage, as shown in Fig. 11. One can see that 
1 / K ab < except for a narrow region whose location 
varies from 1/4 to as v deviates from 1. Therefore, al- 
though the exact value of f p b is not known precisely, we 
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may conclude that the system eventually flows to a CI SO 
phase. 




FIG. 11. Typical behavior of the Luttinger exponent of 8^ 



when 6 P , is frozen. 

ab 



2. Umklapp and "special" processes 



Umklapp scattering and "special" processes which pre- 
serve momentum but exist only for specific values of the 
Fermi wave-vectors 42 are allowed on special lines in pa- 
rameter space, as shown, for example, in the phase dia- 
grams of the 2-leg ladders. In the following we describe 
all the Umklapp and "special" terms in the presence of 
two pairs of Fermi points. For this purpose we define 



(A5) 



An intra-band Umklapp term exists in the case where 
one of the Fermi points pairs resides at kp i = 7r/2, 



H, 



(i) 

Umklapp 



= dx 2irVi 



+ h.c. 



(A6) 



a,f3 



The RG equation for the single-band Cooper-scattering 
c p ti is altered by this Umklapp process. Denoting the ad- 
dition to the RG equation as Sea we find 



(A7) 



where the dot indicates a derivative with respect to the 
length-scale I, An additional RG equation for the Umk- 
lapp coupling itself is needed, 



4r p v 



(A8) 



with the initial value uu = U /&irvi for the systems that 
we have considered. 

When the sum of kp i and kp. at the right moving Fermi 
points equals 7r, three inter-pair Umklapp interactions 
emerge. 



Hulkiapp = ldx^w(vi + Vj) 

a, ft 



K fl]P K iai/3 + h-C 



+ < 



rL iaj{3 Ti -iajf3 

K m K L 

1 ^iaj{3 r ^jaif3 



h.c. 
h.c. 



+ i <-> j 



(A9) 



The same scattering processes also exist, as momentum- 
conserving "special terms" , when the two wave- vectors 
equal each other. These interactions induce additional 
terms in the Cooper- and forward-scattering RG equa- 
tions 



p _ 



16(u?.) 2 + 4(^) 2 + 4(^) 2 +4<u 



16«,) 2 



b 



8u£(u: 
32<.«+<,) , 



(A10) 



where otij = + Vj) /4viVj. The RG equations for the 
new couplings are 



V-ij ^fij u ij "1" c ij{ u ij U ij) ^ C ij( U ij U ij) ' 

u% = u£(4c& + c%) + u*.(c> + 4.) - + c?.) 

+ u ij{^fij + o-^i^ _ ~2 U ij( C ii + C jj) ! 
^ij — u iji~^ c ij c ij) ~ U ijfij u ij( c ii c jj) 

+ \ut j {c° i + c%) + 2uy? j -\u%ff j , (All) 

with the initial values zt°- = U/8n(vi + Vj), — 
U/4n{vi +Vj), ufj = -U/4n{vi + v 3 ). 

More inter-pair Umklapp processes exist when the dif- 
ference between kp i and kp. at the right moving Fermi 
points equals tt, 



H 



(3) 



U mklapp 



y »3 3 l 



(A12) 

The additions to the Cooper- and forward-scattering RG 
equations are in this case 



"u«) 2 + T7T«y«) 2 ' 



1 



P ^2 



-«) 



2<<- 



3 

1 



«) 2 



(A13) 



and the RG equations for these Umklapp couplings are: 
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u ij - u ij( c ii + c j 3 2 fij) 16 u ij( c S + c jj 2 fi]) ' 



(cS+c^.-2/5) + «^i+4- 2 /«) 



2 U ij( C ii C jj + 2 A?) ' 



(A14) 



with the initial values = U/A-K{vi + Vj), ufj — 
-U/n(vi +vj). 

In the diagonal 2-leg ladder, interactions that involve 
an odd number of operators per band are allowed on spe- 
cific lines. For example, when the sum of and 3k Fj at 
the right moving Fermi points equals zero or equals 2n, 
we find: 



H, 



uLlapp = / dx J2 n ( Vi + W J"K 

+ K^K? a]fi + h.c 



(A15) 



This coupling induces the following modifications to the 
RG flow equations 



6c p jj = 2(3 ij (uf j ) 2 , 



(A16) 



where /3jj = («j + Vj)/2vj. The additional RG equation 
is 



l ij ^ u ij (°jj + c ij + 



(A17) 



with the initial value = U /2n(vi + Vj). 

In the case that the difference between 3/c_F; and fcpv at 
the right moving Fermi points equals 2tt, the Hamiltonian 
is: 



i.R 



(A18) 



The additions to the Cooper- and forward-scattering RG 
equations in this case are 



5c p 



<Jc& = 2/3 ji « j ) 2 • 
The RG equation for the new coupling is 

u ij — u ij {^ c ii — fij — fij) j 

with the initial value u|- = U/2n(vi + Vj). 



(A19) 



(A20) 
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